Random Walk and Broad Distributions on 
Fractal Curves 

Seema E Satin A.D.Gangal 
January 18, 2013 

Abstract 

In this paper we analyse random walk on a fractal structure, specifi- 
cally fractal curves, using the recently develped calculus for fractal curves. 
We consider only unbiased random walk on the fractal stucture and find 
out the corresponding probability distribution which is gaussian like in 
nature, but shows deviation from the standard behaviour. Moments are 
calculated in terms of Euclidean distance for a von Koch curve. We also 
analyse Levy distribution on the same fractal structure, where the dimen- 
sion of the fractal curve shows significant contribution to the distrubution 
law by modyfying the nature of moments. The appendix gives a short 
note on Fourier transform on fractal curves. 

1 Introduction 

Random walks in one dimesion is a well established topic [1] [2 El H] , with several 
applications in physics including transport phenomena and diffusive behaviour. 
The discrete simple random walk in 1-dimension is described by a particle taking 
steps, each towards the right or left with a specified probability. For unbiased 
walks the probability is same for right or left step 3^ . In the continuum limit the 
number of steps N goes to oo, the step size 5 goes to zero, such that NS remains 
fixed. The random variable is rescaled to give the probability distribution [21 [3] , 
showing the diffusive behaviour. 

While studying transport properties on disordered systems, one encounters 
subdiffusive behaviour [5 . Some physical examples are, NMR diffusometry on 
percolation stuctures ^14 , motion of a bead in polymer network 15 etc. In 
general the diffusion process can be expected to depend on both the geometry 
of the medium and the process. While various aspects of dependence on pro- 
cesses (e.g. Gaussian, Levy etc) and ordinary geometries are well studied, the 
dependence on fractal geometries of media is relatively unexplored. 

It is expected that on fractals, subdiffusion can be due to the geometry of the 
underlying fractal stucture on which the process takes place. In this chapter, 
various well established approaches in random walk problems are extended to 
fractal geometry. We study such a behaviour by considering a random walk on a 
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fractal curve . In particular, we present the analytical and numerical results for 
walks on a von Koch curve in R^. Random walk on such a curve is performed 
by taking a step either forward or backward along the curve itself with equal 
probability if the walk is unbiased. Here we consider unbiased walk only. 

Such a random walk in the continuum limit leads to study of processes like 
diffusion on the Fractal curve F. Hence one can study the diffusive behaviour of 
particles, performing random motion on the fractal curve, by using this formu- 
lation. Appropriate and exact expressions can be obtianed, for a class of walks, 
to describe physical processes by using this model. The simplest case is that of 
obtaining a probability distribution, which is same as the solution obtained by 
solving a diffusion equation on a fractal curve as in [T7] . The study of a discrete 
random walk on a fractal curve F is another approach towards studying such 
transport properties on these kind of structures. 




Figure 1: Example of a fractal curve. 9 and 9' denote any two points on the 
curve, and I gives the Euclidean distance between these two points 

There are several approaches to study anomalous diffusion. Some commonly 
used frameworks are fractional Brownian motion [B], the continuous time ran- 
dom walk |/7j , fractional diffusion equations [51 [S] , generalized Langevin equation 
and Fokker-Planck equations [TUJ [TTl [TH [13] etc. In most of these approaches 
the assumptions for validity of simplest form of Central limit theorem are not 
satisfied _5j. 

In the present case, the Central limit theorem is not violated. The departure 
from Gaussian distribution or, the normal diffusive behaviour, is explicitly due 
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to the underlying fractal geometry of space and not due to correlations or long 
tailed statistics. 

Here, we consider a random walk on a fractal curve(von-Kocli like), in R". 
The calculus on such a curve is developed in [T^ , where the curve F has fractal 
dimension a. 6 denotes a point on the curve and J (9) = Sp{u) denoted the 
mass function, which is the mass of the fractal covered upto point 6, u is the 
parameter of the curve. A short review of the Calculus is given in Appendix 2. 

Random walk on such a curve is performed by taking a step either forward 
or backward along the curve itself with equal probability if the walk is unbiased. 
In this paper we consider unbiased walk only. 

In various sections of this paper we find analogues of simplest aspects of 
random walk problems in 1-D extended to fractal curves. 

2 Discrete Random Walk on a fractal curve 

We consider the special class of random walks on a fractal curve where the 
walker covers a fixed mass A in each step. If the fractal curve, F reduces to 
a real line, then the random walk on it reduces to an ordinary simple random 
walk as analysed in [18]. 

Let C{N, 9, 9') be the number of walks that start at a point 9 on the fractal 
curve F with 7-dimension, a ( note that a > 1 in general, particularly here 
1 < a < 2 ), and end at another point 9' on the same fractal curve, at the N'^^ 
step. Then 

C(7V; 61, 9') = C{N - 1; 6*, .r^{J{9') - A)) + C{N - 1; 9, .r^{J{9') + A)) (1) 

Equation ([T]) is the recursion relation for a random walk on F. 

When N i.e number of steps is large, and A small, equation ([T]) can be 
approximated by a differential equation involving the i^"-derivative. This is 
carried out in the following. 



C{N + 1; 9, 9') = {C{N- 9, J-\J{9') - A)) + C{N; 9, J-\J{9') + A)) - 2C{N- 9, 9')} 
+2C{N,9,9') 

(2) 

Taylor expanding (as given in appendix) first two terms around J{9') in powers 
of A, we get 

C{N + 1; 9, 9') = {A^{D'^^, fC{N; 9, 9') + ^ A^p^^, fC[N- 9, 9') + ...) + 2C{N; 9, 9') 

K. A^{D'^^,fC{N;9,9') + 0{A'^) + 2C{N;9,9') (3) 

(4) 

C{N + 1; 9, 9') - 2C\N, 9, 9') w A^iD^^p , fC{N- 9, 9') (5) 
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Each step has equal probabihty in forward or reverse direction on the curve. 
The total number of steps being N, there are 2^ ways of performing an N step 
random walk. Thus, P{N,9,9'), the probability of a such an TV-step random 
walk, starting at 9 and ending at 9' is given by 

P{N,9,9')^^C{N,9,9') 

Hence we can replace C{N, 9, 9') by 2^P{N, 9, 9') and C(iV+l, 6*, 9') by 2^+^P{N- 
1,9,9'). 

Thus, equation ^ leads to 

P{N + 1, 9, 9') - P{N- 9, 9') « ^{D'^F,, fP{N; 9, 9') 

We now assume N to be large and P{N,9,9') to be slowly varying function of 
N. Thus we can write the above equation in the form 

A.P(^N;9,9') = ^iDF„rPiN,9,9') (6) 

The solution of the above equation can easily be obtained by using conjugacy 
between the _F"-derivative and ordinary derivative as in [16 . This gives. 

Now, one can go to the continuum case, where one can set t = Nt, N being 
the number of steps in the discrete walk and t being the duration between 
the consequtive steps. Thus the expression for probability distribution in the 
continumm case, which is of the form given by equation (1111) in the later section 
can be obtained and used readily. 



3 Moments 

We now consider the Euclidean distance L{9) = L{w{u)) = |w(u)| upto a point 
9 — w(u) from the origin on the fractal curve F. 

Consider the probability distribution which behaves as 

P{9) = P(w(m) =9)^ exp(-S'^(u)2) where S'^(m) = J{9) (8) 

with appropriate normalization. Using this form for the probability distribution 
we calculate the first two absolute moments on the fractal curve as follows: 

<L>= f L{9)P{9)d%9 (9) 

J C( — oo.oo) 

and 

< >= / L^{9)P{9)d%9 (10) 

J C( — oo.oo) 
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We evaluate the above two moments in terms of Euclidean distance, using 
the explicit form 

as obtained in chapter [17] , which is the continuum case of equation[7|as N ^ oo 
and A -J> (TV A being held fixed). 

Heuristic Calculation of Absolute Moments 

Now, for the given gaussian like probability, eauation (|lip . we calculate the 
moments in terms of Euclidean distance, for random motion on a von Koch 
curve. From the graph of Rise function and Euclidean distance, for a von Koch 
curve as given in [16], it is clear, that ci[L(w(u))]" < Sp{u) < C2[L(w(u))]". 
Thus we can approximate Sp(u) by L", i.e 

S^{u) ^ [L(w{u)r 

or in short hand notation 

Thus, the first moment can be calculated analytically by using an intuitive 
replacement for the F" -integrals in equations (|9]) and (fTOl) . Sp is thus replaced 
by L" and dp9 is replaced by dL°', keeping in view the way _F"-integrals are 
defined in [THj . 

Thus, 

On substituting /\/2At ~ z we reduce the integral to the fohowing 

poo 

< i >= Const.ii/2« / z^/°'exp{-z^)dz (13) 
Jo 

which implies 

similarly we can obtain the second moment in terms of the Euclidean distance 

/ ^^exp(-(|^ 
leading to the behaviour 

< >- t^/°' (16) 

In the case of von Koch curve 1/a — 0.792 and l/2a ~ 0.396. These 
expressions are in accordance with the numerical results presented in fig ([2]) and 



< 2-i= / L'eM~^-^)dL" (15) 
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Figure 2: Plot of log < L > (Y-axis) vs. log{t) (X-axis). The cross denotes 
the value (calculated by performing F"'- integration numerically), to which the 
straight line is fitted with slope 0.403 , to extract the value of exponent in 
equation ([T4|. for a von Koch curve. 

([3]). In these figures a plot between moments and time is shown. These plots 
have been obtained by performing the integrals in equations ^ and (|10p 
numerically, for the von-Koch curve. From the plots one obtains 

and 

Here we see a beautiful match between the expected values and numerical 
calculations within numerical accuracy. The plots also comfirm the expected 
behaviour of Sp{u). 

4 Broad Distributions 

We have derived in section ^ , the expression for probability distribution for a 
discrete, simple random walk on fractal curve. Now we examine random walks 
with some different (non Gaussian) probability distributions for the individual 
steps. 
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Figure 3: Plot of log < > (Y-axis) vs. log{t) (X-axis). The cross denotes 
value (calculated by performing .F^-integration numerically) to which straight 
line is fitted with slope 0.802 to extract the value of the exponent in equation 
([16] ), for a von Koch curve. 

Review of Broad Distributions 

We begin by summarizing some relevant results for Broad distributions from 



For a 1-D random walk in ordinary space, let X^q denote the sum of inde- 
pendent random variables Z„, i-e 



where In is the n*^ step size. Let us denote p{l), the probability distribution for 
^ (individual steps in the random walk) which is broad, that is ,it decreases for 
large I as with /i > 0. We consider two cases, 

• For < /i < 1 , Xiq , the first moment (l) is infinite. 

• For 1 < /i < 2 , the first moment {1} is finite and the second moment {P) 
is still infinite. 

It is well known that limit distributions of the sum are defined by their 
characteristic functions which are given by Lp_f^{k) in Fourier space, i.e k- space. 



El- 



N 




(17) 
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when fi and /3 are parameters, /3 being the degree of asymmetry. We consider 
the case /3 = for the stable distribution io,/i- In this case large positive and 
negative values of In occur with equal frequency. Then 

LoA'^) = e""'" (18) 

Its Fourier transform being 

We set Z ^ Zn where,ZAr Xn/N^'^" or [Xn- < I > N)/N^'^' as iV ^ oo 
,for 0</i<lorl</i<2 respectively. All the important results for such a 
case are given in |5j. 



Broad Distributions on fractal curves 

Now we consider an analogous approach for the random walk on fractal curve. 
In particular we take parametrizable fractal curve, the von-Koch curve with 7- 
dimension equal to logA/logS as described in [TB] . 

We consider the following probability distribution in the limiting case, when 
N ^ 00. We now abuse the above notation for obvious reasons to denote 
Lo,ij,{ip) as the distribution in Fourier space of fractal curve F, where ip = w(fc) 
and Sp{k) — J{ip) as described in appendix. Then 

Lo,M=eM-\Jmn (20) 

In what follows we will use L to refer to quantity on the fractal curve (space) 
and L to refer to the quantity in the real line. 

The inverse Fourier transform (as defined in the appendix), for the above 
distribution is given by 

LoA(^) ! exp(*J(^)J(0) - I J(^)r)d^V (21) 



27r 



C( — 00,00) 



The above integral in view of conjugacy between F°' integral and Riemann 
integral becomes: 

UAV = J{Q)) = ^ / exp(zfcy - \k\^)dk (22) 
For n — 2 we obtain the gaussian distribution, case: 

Zo,^(y = J(0)) = -^exp(-^) (23) 

Applying the fractalizing operator to the above equation we get back the ex- 
pression for distribution on fractal curve. 

i„,,(.,.-^o.p(-MMdl!, (24, 
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Similarly for the Cauchy case when fi — 1 we get 

1 



Lo,fi — 



7ril + {J{e)^S^{u)}^) 



Now, one can then discrctize the integral in equation (|22p by replacing dk in 
terms of equal steps of size Afc, as follows: 

oo 



where km — mAk. 

For values of /i other than 1 and 2 , the expansion for large arguments is 
given by [S] 

^-^ ml 

m— 1 

The leading term of which is 

L^^oiy = J{0)) = {^)-^V-^^+^^T{l + ^i) s\n{^:^i/2) 
Applying the fractalizing transformation to the above we get the result: 

IMS) = (7r)-i5^(w)-(^+i)r(l + /i)sin(^/i/2) (25) 
- (26) 

for large values of J{0) = Sp{u) or L where L denotes the Euclidean distance. 

One can see that < L > is infinite for /x < 1/a and is finite for /i, > 1/a for 
these large values of L. Also the second moment < > diverges for /i < 2/a 
and is finite for /i > 2/a for large values of L. Hence we see a scaling down of 
the absolute moments for Levy distrubutions due to the fractal structure of the 
underlying space. 



5 First Passage Time 

The first passage time is a standard problem in Statistical Physics [121 [3] . Here 
we extend it to first passage time on fractal curves. 

Random walk on a fractal curve can be described by forward and backward 
steps along the curve that cover equal mass on the fractal curve as in section ^ , 
so that there is one- to one correspondence between random walk on a straight 
line (the parameter u) and that on the curve F. 

The first passage time, or the time required to reach a certain point 9 = w(u) 
on the curve F, for the first time can be defined as follows: 
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Definition 1 The first passage time T{9) to a point 9 € F is given by 



T{e) = inf{i : e(t) ^ 9} (27) 

where Q{t) is the position of a random walker on F at time t, which started at 
a point 6{) at time i = . 

Now we explore the relation between this first passage time and the mass covered 
by the random walker on the fractal curve in this time. 

When the walker performs forward and backward steps along the curve, let 
the mass covered in the n}^ step be given by 7|^(F, u„, m„_i), corresponding 
to parameter value lying between u„_i and u„. The resultant displacement 
( considered on the 5^-axis), of the particle performing a random walk of M 
steps, some of which may be in forward direction and rest in the backward 
direction on the curve, so as to reach an arbitrary point w(ufc) on F , is given 
by sum of over M steps. For k < M such that M = k + r, and step size A, 
for every sequence {ui, U2, . . . , wm}, such that, 

SUum) = kA (28) 

Let the time duration between consecutive steps on the fractal curve F, be 
denoted by t. We calculate the maximum mass that a random walker covers in 
a given time t, such that t = Mr, where M is the total number of steps taken 
in time t. 

Since k = M — r wc can write equation (|28p in the following form: 



5^(uM) = (--r)A (29) 

T 

It can be clearly seen that Sp{uM) is maximum when r = on the rhs of the 
above equation. Hence 

SUuM)\ma^ = -A (30) 
T 

Next we calculate the minimal time required by a random walker to cover 
certain mass on the curve. 

Conversely, let t — Mr be the time taken for M steps on the curve -F, also 
let M = fc + r as given above, then k ~ -^r^, and we can write equation [28] as 

/ — TT 

SUum) A (31) 

r 

Thus, 

TSf{uM) 

t =^ rr ^ 

A 

Hence we see that t is minimum when r — 0, Hence 

tmin = -^Sp{um) (32) 
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Figure 4: Plot of Lmax (Y-axis) vs tmin (X-axis) for random walk on the von 
Koch curve, obtained numerically. 

An interesting quantity, which can be analysed is the maximum Euclidean 
distance on a fractal curve, covered by the random walker in a minimum time 
t. In general these will depend on the geometry of the particular fractal curve 
of interest. The relation between Euclidean distance on a fractal curve and the 
staircase function [16^ can only be obtained numerically, even for a simple fractal 
curve like the von-Koch curve. A plot between t„iin minimum time required to 
reach a point and Lmax (maximum Euclidean distance) calculated numerically 
for a von Koch curve, is shown in fig One may note a step- like behaviour 
in this figure. We comment on this in the results. 

6 Results 

It is clear from the figures ^ and ([3]) that a subdiffusive behaviour is addressed 
by the gaussian-like random walks on the von Koch curves. 

Eq. (|25p shows the behaviour of the stable law o for values of 9 for 
which Sp(u) is large. This behaviour differs from the ordinary law for large 
deviation, when the underlying space is not a fractal and when the stable law 
decreases as y^'-^^'-'\ In the case of an underlying fractal space, the exponent 
gets multiplied by a factor of a the dimension of the fractal curve itself. Thus we 
see that from a Euclidean perspective, the behaviour gets a direct contribution 
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from the exponent of the distribution and the dimension of the curves. This is a 
striking difference. When taking straight hne approximations on fractal curves, 
the dimension of the fractal curve contributes only indirectly. While using the 
above method, the fractal dimension plays a direct role in changing the nature 
of heavy tails of the distribution. 

In fig. (jH) we see a step-like behaviour, this explicitly shows that the Eu- 
clidean distance on the curve remains constant for a short while and then rises. 
There is a largest period of constant Euclidean distance in the figure, this de- 
pends on the geometry and the initial position of the random walker on the von 
Koch curve. Here we calculate the Euclidean distance from the origin or initial 
point (0,0) on the curve in . 



Appendix 1 

The Fourier Transform 

From the definition of conjugacy |16] 

miSf{u))^fiwiu)) (33) 
The Fourier Transform on the real line, for a function g{v), is defined by 



9iv) = / g{y) ex.p{-ivy)dy (34) 

J —oo 

and the inverse Fourier Transform is 

1 f°° 

9{y) = TT / 9{v) exp{ivy)dy (35) 

In the case of a parametrizable fractal curve F, which is obtained by a fractal- 
izing transformation on an interval of the real line, we propose that the Fourier 
space can also be obtained by the same fractalizing transformation on the in- 
terval of a real line. The interval may be (— oo, oo). 

Let g = (j)[f] and g = (j)[f] also v = Sp{k) and y = Sp{u). 

We use the notation J{0) = Sp{u) and J('0) = Sf{k), where 9 = w{u) and 
tjj = w(fc). 

Then taking Fourier transform of the LHS of equation pSI) one can write 



(/.[/](!; = J(^)) = / ^[f]{y = J{9))exp{-2yv)dy 

J —OO 

= / f{9)exp{-iJ{9)v)d'^e (36) 

Jc(-oo.oo) 



where 



C(— oo,oo)= lim Cia,b) 
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and 

f{e)exp{^iJ{e)v)d'^9 (37) 

(7( — 00,00) 

Comparing equations and (f57|) . we can write 

Also one can define the action of in Fourier space as 

^[f]{v = JW) = m (38) 
Now using conjugacy one can rewrite equation p6p as 

fW^ f /(^^)exp(-^J(0)J(V))d^0 (39) 

C( — 00,00) 

Similarly, inverse transform of the above can be obtained from equation psp . 
which can be written as 



f{0) = ^l f{ij)exp{tJie)j{i:))d'^e (40) 

^TT ./c(-oo,oo) 



Appendix 2 

Review of Calculus on Fractal Curves 

For a set F and a subdivision -P[a,6], a < b, [a, 6] C [oq, 6o] let w : [a, b] — > 
then we define the mass function as follows: 

7"(^^,a,6) = lim inf y Ij^^lhlll^J^ (41) 

where | • | denotes the euclidean norm on R", 1 < a < n and P[a,b] = {a = 
to, ■ ■ ■ ,tn = b}. 

The staircase function, which gives the mass of the curve upto a certain 
point on the fractal curve F is defined as 

[ -7"(i^,M,po) u<po 

where u € [oq, 60] • 

A point on the curve w(m) = 9 and Sp{u) = J{0). The derivative is 
defined as: 
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The F"-integral is also defined and is denoted by 

/ fiO)d^O (44) 

JC{a,b) 

where C{a, b) is the section of the curve lying between points w(a) and w(6) on 
the fractal curve F. • 
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